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According to the theoretical results, the quantum searching al- 
gorithm can be generalized by replacing the Walsh-Hadamard(W-H) 
transform by almost any quantum mechanical operation. We have im- 
plemented the generalized algorithm using nuclear magnetic resonance 
techniques with a solution of chloroform molecules. Experimental results 
show the good agreement between theory and experiment. 
PACS number(s):03.67,89.70 



The quantum searching algorithm was first proposed by Grover [1]. It 
can speed up some search applications over unsorted data and has been 
become a hot topic in quantum information [2~5]. The surprising result 
of the algorithm is that searching a particular item in unsorted list of N 
elements only requires 0{\/N) attempts [2]. The algorithm needs the prop- 
erties of Walsh-Hadamard(W-H) transform and has been implemented using 
NMR techniques [6]. Grover generalized the algorithm in theory, in which 
the Walsh-Hadmamard transform can be replaced by almost any quantum 
operation [7]. 

The generalized algorithm can be posed as follows. If a unitary operator 
U is applied to a quantum system in an initial basis state I7 >, the system 
lies in a superposition J^i^ali where aj =< z|f/|7 >= Ui-y. The amplitude 
of reaching the target state |r > is Ur-y, and the probability of getting the 
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system in state |r > is |C/r7l ^ measurement is made. Therefore, it will take 
at least 0(1/|[/t-^P) repetitions of U and the measurement to get state |r > 
for one time. In contrary, according to the generalized searching algorithm, 
in 0{l/\Ur^\) applications, operator Q transforms I7 > into t/~^|r >. After 
U is apphed, the system lies in state |t >. Q is defined as Q = —I^U~^IrU , 
where /j = I — 2\i >< i|, and I denotes unit matrix. If 7j is applied to a 
superposition of states, it only inverts the amplitude in state \i >, and leaves 
the other states unaltered. It has been showed in theory that the number of 
repetitions of Q required to transform I7 > into |r > is 7r/4 16^1-71 if I^/ttI 1. 

It can be proved that for a two qubit system, the number of repetitions 
of Q is also 0{l/\UT'y\) without the condition that \Ut^\ <C 1. Particularly, if 
\Urr^\ = 1/2, the number of repetitions is 1, and the probability of getting the 
target state is 1. In this paper, we will implement the generalized quantum 
algorithm in this case on a two qubit NMR quantum computer [8] . 

Our experiments use a sample of Carbon- 13 labelled chloroform (Cam- 
bridge Isotopes) dissolved in d6-acetone. Data are taken at controlled tem- 
perature (22''C) with a Bruker DRX 500 MHz spectrometer (Beijing Nor- 
mal University). The resonance frequencies vi = 125.76 MHz for ^^C, and 
1^2 = 500.13 MHz for ^H. The coupling constant J is measured to be 215 
Hz. If the magnetic field is along 5-axis, and let ^ = 1, the Hamitonian of 
this system is 

H = -2nuill - 2t:v2II + 'inJilll, (1) 

where I^{k — 1, 2) are the matrices for i-component of the angular momen- 
tum of the spins [9]. In the rotating frame of spin k, the evolution caused by 
a radio- frequency(rf) pulse along x or y-axis is denoted as Xk{'^k) = e^'^''^^ 
or Yi:{—ipk) = e~^''^''^y , where ipk = Bijktp- -Bi, Jk and tp represent the 
strength of magnetic field, gyromagnetic ratio of spin k and the width of the 
rf pulse, respectively. The pulse used above is denoted as [f/?]^ or [— V'ly- The 
coupled-spin evolution is denoted as 

[t] = e-^^^-'Ji^l^ (2) 
where t is evolution time. The pseudo-pure initial state 
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is prepared by using spatial averaging [10], where | |>j. denotes the state 
of spin k. For convenience, the notion | |'>i | |'>2 is merged into | |"f> 
and the subscripts 1 and 2 are omitted. The basis states are arrayed as 
I TT>)| Ti>)| iT>) I ii>- III matrix notion, /j is a diagonal matrix with 
all diagonal terms equal to 1 except the ii term which is -1. For example, if 
\i >— I ii>, Ii is expressed as 
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U is chosen as three forms 



Ui = 1/2 
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and 



U2 = 1/2 
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where C/i is the W-H transform used in reference [6] (up to an phase factor). 
It is easily proved that UQ can transform state I7 > into state |t >. The 
irrelevant overall phase factors are ignored. 

The evolution of the system can be represented by its deviation density 
matrix pa[11]- For the heteronuclear system, the following rf and gradient 
pulse sequence [a\l - [grad], - [7r/4]^ - 1/4 J - - - 1/4 J - [-7r/4]], - 
[gra(J\z transforms the system from the equilibrium state 



PAeq = 111] + I2I. 
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to the initial state required in experiments 

PAo = -(4V2 + ll/2 + Illl) = -1/4 
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where a — arccos(7i/272), and [grad]z denotes gradient pulse along z-axis. 
The pulses are applied from left to right. The symbol 1/4 J denotes the evo- 
lution caused by the Hamitonian H for 1/4 J without pulses. The initial state 
can be used as the pseudo-pure state | 'f'|> in NMR quantum computation 
[12]. By applying pulse [vr]^, [vr]^, or [7r]J^^[7r]^, the pseudo- pure states corre- 
sponding to states I tJ,>, I |t> or | ll> can be gotten. They are written 
as 

PAi = -{11/2 - 11/2 - Illl) (10) 

PA2 = -{-11/2 + 11/2 - Illl) (11) 

p^^ = -{-ll/2-ll/2 + llll) (12) 

We also represent the four pseudo- pure states above as | tt>,| Ti>)| iT>) 
and I ii>, respectively. According to reference [6] , we find that 

7o = yi(7r/2)F2(7r/2)Xi(-7r/2)X2(-7r/2)ri(-7r/2)F2(-7r/2)[l/2J], 
h = Fi(vr/2)F2(vr/2)Xi(7r/2)X2(-7r/2)ri(-7r/2)r2(-7r/2)[l/2J], 
h = Fi(7r/2)F2(7r/2)Xi(-7r/2)X2(7r/2)ri(-7r/2)r2(-7r/2)[l/2J] and 
h = Fi(7r/2)F2(7r/2)Xi(7r/2)X2(7r/2)Fi(-7r/2)F2(-7r/2)[l/2J]. The time 
order is from right to left. It is easy to prove that 
U, = Xi(7r)X2(7r)Fi(-7r/2)y2(-7r/2), 
U2 = Fi(7r/2)r2(7r/2), and 
f/3 = Xi(7r/2)X2(7r/2). 

In our experiments, each NMR spectrum of spin k is obtained by a spin- 
selective readout pulse [— 7r/2]^. The relative phases of signals are meaningful 
because all experiments are acquired in an identical fashion [13]. At first, we 
prepare the four pseudo- pure states described in equations(9)-(12). For the 
system in a pseudo- pure state, only one NMR peak appears in a spectrum 
if a selective readout pulse is applied. Fig.l shows the experimental MNR 
spectra when the system lies in various pseudo- pure states. Figs. la, b, c and 
d are spectra when the system lies in pseudo- pure states | tT>)| Ti>)| iT> 
and I ii>, respectively. In each figure, the main spectrum represents the 
spectrum of ^^C and the spectrum is a smaller inset. For each experiment. 
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the initial state I7 > is pseudo- pure state | tT>- U is selected as Ui, U2 or C/3. 
The searching state |r > can be any pseudo- pure state. UQ transforms state 
I7 > into state |r >. Therefore, the number of application of Q is 1. The 
searching results on condition that U = U3 and |r >= | tt>, | ^i>, \ iT>, 
or I J,J,> are shown in Fig. 2. Comparing Figs. 2 a, b, c and d with Figs. la, b, 
c and d respectively, we confirm that the system truly lies in the target state. 
In order to illustrate that the searching results for f/3 are the same as Ui or 
U2, Fig. 3 shows the searching results on condition that |r >= | ti>) ^-iid 
U = Ui (shown by Fig. 3a) or U2 (shown by Fig.3b). It can be found from 
various spectra that the experimental errors are not lager than 5 percent 
expect only two spectra. The errors mainly result from the imperfection 
of pulses, inhomogeneity of magnetic field and effect of decoherence. We find 
that if each [7r/2]^ (</'=x or y) in Jj is replaced by [— 7r/2]^, and each [— 7r/2]^ 
is replaced by [7r/2]^, the results remain the same. This fact can be used to 
simplify pulse sequences and reduce experimental errors. 

In conclusion, we demonstrate the generalized quantum searching algo- 
rithm by replacing W-H transform by different transformations. Because the 
number of repetition of operator Q is determined by the element \Ur-y\ and 
these three U transformations have the same lUr-yl, it is not surprised that 
the numbers of repetitions of Q are all 1. Compared with temporal labelling, 
the spatial averaging used to prepare initial state shortens experiment time. 
It only takes about 2 minutes to finish one experiment. The searching results 
can be directly read out from spectra, so that the steps of recording areas of 
peaks are avoided. These facts simplify the process of experiments and make 
the searching results easy to observe. 

This work was partly supported by the National Nature Science Foun- 
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Figure Captions 



1. The NMR spectra of ^^C(main figures) and (smaller insets), when 
the two-spin system lies in various pseudo- pure states via readout 
pulses selective for ^^C (main figures) and ^if (smaller insets). The 

amplitude has arbitrary units. When the system lies in a pseudo-pure 
state, only one NMR peak appears in the spectrum of spin k if a readout 
pulse [—7i/2]y is applied. Figs, la, b, c and d are spectra corresponding 
to pseudo- pure states | tT>) I ti>; I it>8'iid | ii>, respectively. 

2. Spectra of ^^'C (main figures) and (smaller insets) after completion 
of the generalized searching algorithm and readout pulses selective for 
^^C (main figures) and (smaller insets). U is chosen as C/3, and the 
target state |t > is | tT>) I Ti>) I iT> or I ii>- The corresponding 
searching results arc shown by Figs. 2a, b, c or d. By comparing Figs. 2a, 
b, c and d with Figs. la, b, c and d respectively, we confirm that the 
system is truly in the target state. 

3. Spectra of ^^C (main figures) and (smaller insets) after completion 
of the generalized searching algorithm and readout pulses selective for 
^^C (main figures) and (smaller insets) on condition that |r >= 
I |j>, and U — Ui shown by Fig. 3a or U = U2 shown by Fig. 3b. 

[Figure 1 about here.] 

[Figure 2 about here.] 

[Figure 3 about here.] 
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